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Abstract 


Organometal halide perovskites are promising solar-cell materials for next-generation pho¬ 
tovoltaic applications. The long carrier lifetime and diffusion length of these materials make 
them very attractive for use in light absorbers and carrier transporters. While these aspects 
of organometal halide perovskites have attracted the most attention, the consequences of the 
Rashba effect, driven by strong spin-orbit coupling, on the photovoltaic properties of these 
materials are largely unexplored. In this work, taking the electronic structure of methylam- 
monium lead iodide as an example, we propose an intrinsic mechanism for enhanced carrier 
lifetime in 3D Rashba materials. Based on first-principles calculations and a Rashba spin-orbit 
model, we demonstrate that the recombination rate is reduced due to the spin-forbidden tran¬ 
sition. These results are important for understanding the fundamental physics of organometal 
halide perovskites and for optimizing and designing the materials with better performance. 

The proposed mechanism including spin degrees of freedom offers a new paradigm of using 
3D Rashba materials for photovoltaic applications. 

The organometal halide perovskites (OMHPs) have attracted significant attention due to the 
rapid increase in their photovoltaic power conversion efficiency. In the past 2 years, the reported 
efficiency of OMHP-based solar cells has almost doubled from 9.7%®to over 20%J 2 4 making 
OHMPs very promising for low-cost and high-efficiency photovoltaics. Methylammonium lead 
iodide, CH 3 NH 3 PM 3 (MAPM 3 ), and other closely-related hybrid perovskites such as Cl-doped 
and Br-doped MAPM 3 (MAPbl 3 _ x Cl x and MAPbl 3 _ x Br c ). (NtbbCHPb^ (formamidinium lead 
iodide, FAPM 3 ), and Sn-doped MAPM 3 (MAPb x Sni_ x l 3 ), all display band gaps (1.1 to 2.1 eV) 
in the visible light region, favorable for photovoltaic applications . 5 13 The class of materials also 
possesses strong light absorption, fast charge generation and high carrier mobility. 14 15 In partic¬ 
ular, exceptionally long carrier lifetime and diffusion length have been observed in MAPM 3 and 
MAPbl 3 _ x Cl x , making them better solar cell candidates than other semiconductors with similar 
band gaps and absorption coefficients . 16 18 

Intense research has been directed toward understanding and further enhancing the long car¬ 
rier lifetime and diffusion length in OMHPs. Previous studies reported a relatively low defect 
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concentration in MAPbl 3 , 19 ' 23 which reduces the scattering centers for nonradiative charge car¬ 
rier recombination. Recently, it has been suggested that the spatial carrier segregation caused 
by disorder-induced localization 24 or domains acting as internal p-n junctions 25 27 may reduce 
the recombination rate. Though many of the OMHPs are 3D Rashba materials driven by strong 
spin-orbit coupling (SOC) and bulk ferroelectricity, 28 31 the effects of spin and orbital degrees of 
freedom on photovoltaic applications are largely unexplored beyond band gap engineering. 29 In 
this work, we focus on an intrinsic mechanism for the enhancement of long carrier lifetime due 
to the Rashba splitting. Using first-principles calculations and effective models, we find that the 
Rashba splitting arising from SOC under inversion symmetry breaking can result in spin-allowed 
and spin-forbidden recombination channels. The spin-forbidden recombination path has a sig¬ 
nificantly slower transition rate due to the mismatch of spin and momentum. The spin-allowed 
recombination path, though kinetically favorable, can be suppressed under appropriate spin tex¬ 
ture due to the low population of free carriers. Taking the electronic structures of MAPbl 3 under 
various distortions as examples, we show that the proposed mechanism is possible under room 
temperature, and is potentially responsible for the long carrier lifetime in OMHPs. This spin- 
dependent recombination mechanism highlights the possibility of using 3D Rashba materials for 
efficient photovoltaic applications. 

Fig. |T| illustrates the mechanism for enhancing the carrier lifetime in a generic 3D Rashba 
material. The strong spin-orbit coupling effect from heavy elements (e.g., Pb, Sn, I and Br) and 
the polar distortion (e.g., aligned molecular dipoles in OMHPs) give rise to the Rashba effect, 
which lifts the two-fold degeneracy of bands near the band gap. Near the band gap, the spin 
degeneracies of the conduction and valence bands are lifted, giving rise to “inner" and “outer" 
bands with opposite spin textures, characterizing spin rotation direction as “clockwise" (x = — 1) 
and “counterclockwise" (X — +1) (Fig- []}• The photo excitation process creates free electrons 
and holes, which can quickly relax to band extrema in the presence of inelastic phonon scattering. 
When the spin textures of conduction band minimum (CBM) and valence band maximum (VBM) 
are opposite, the radiative recombination of C^=_i —* V^. =+ i is a spin-forbidden transition due 
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to the mismatch of spin states. This prevents rapid recombination as the photon-induced spin-flip 
is a slow process . 32 Moreover, the minimum of the C^=_i band is slightly shifted compared to 
the maximum of the V^ =+ i band (momentum mismatch). This creates an indirect band gap for 
recombination, which further slows down the recombination process due to the requirement of a 
phonon with the right momentum . 33 In the following, we use the terms favorable and unfavorable 
relative spin helicity to describe cases when the VBM and CBM have opposite and aligned spins, 
respectively. 

The spin texture and carrier population of the CBM and VBM play key roles in enhancing 
carrier lifetime. Our first-principles density functional theory (DFT) calculations support the real¬ 
ization of this mechanism in OMHPs. Taking the pseudocubic phase MAPM 3 as an example, we 
explore the carrier dynamics after the photoexcitation. Using Fermi’s golden rule, we calculate the 
inelastic phonon scattering rate (see Methods). As shown in Fig. [ 2 J the relaxation rate (« 10 15 s _1 ) 
increases as a function of carrier energy with sharp jumps corresponding to thresholds for optical 
phonon emission. We find that the organic molecule plays an important role for carrier relaxation. 
The modes that are most important for carrier relaxation are MA translation, CH/NH twisting, and 
CH/NH stretching, and they contribute equally to the relaxation of both electrons and holes. 

Our calculations reveal that the phonon-induced carrier relaxation rate value is many orders 
of magnitude faster than the electron-hole recombination rate (« 10 9 s - 1 )P 4 35 Therefore, the 
carriers will rapidly thermalize and from a quasi-static equilibrium distribution near the CBM and 
VBM. In the ideal case of low temperatures and large Rashba splitting, nearly all free carriers are 
located at the band extrema, and the effects of spin and momentum mismatch on the enhancement 
of the carrier lifetime will be the greatest. This effect is less strong at finite temperatures and 
small Rashba splitting because of the thermal occupation of C^ = i and V^ = _i bands, which opens 
spin-allowed recombination paths such as C 7 - _i —> V^ = _i. In the case of favorable relative spin 
helicity, we investigate this temperature effect by examining the Rashba splitting using a Rashba 
Hamiltonian : 36 
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(2) 


For the parameters m and vr representing the band mass and Rashba interaction respectively, we 
use values obtained from the DFT band structure of fully-relaxed MAPM 3 . We find that the con¬ 
duction band Rashba splitting (0.108 eV) is much larger than the thermal energy scale, while the 
valence band Rashba splitting (0.016 eV) is comparable to the thermal energy scale. Since elec¬ 
tronic correlations are not fully captured in DFT , 13 these values are likely lower bounds of the true 
splitting. These relatively large Rashba splittings are likely to give rise to a significant enhancement 
in carrier lifetime even at room temperature. 

There is an additional factor, arising from the unique features of Rashba band structure, which 
promotes occupation of the band extrema. In contrast to the band extrema of ordinary parabolic 
bands in semiconductors which are points in momentum space, those of Rashba materials are 
one-dimensional rings (Fig. [I]). This leads to an increase in the density of states at low energies, 
resulting in a population of carriers heavily skewed towards the band extrema (Fig. [3^ and b) and 
consequently reduces the overall recombination rate due to the reasons of spin and momentum 
mismatch as discussed above. 

The magnitude of Rashba splitting depends on the amount of polar distortion and the strength 
of the SOC, both of which can be captured by the Rashba velocity parameter vr. We calculate the 
averaged recombination rate 

= y 

he' Jd 3 kn e x (k)n h x ,(k ) 

where = B xx m e x (k)n h x ,{k) is the band- and momentum-resolved recombination rate. The 

spin-mismatch effect is captured by the rate constant B x x i, which is larger when X and x! h ave 
parallel spins than otherwise. Enhancement of density of states enters via the temperature depen- 
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dent electron and hole occupation numbers n e (k) and n h (k), which tend to peak at different k points 
because of momentum mismatch. In order to quantify the effect of Rashba splitting on the recom¬ 
bination rate, we define the unitless lifetime enhancement factor as the ratio (t 1 )v w =o/{t~ 1 )y A> , 
where (t ') v , w is the average lifetime when SOC is taken into account, and (t l ) yff =o refers to a 
calculation where SOC is explicitly set to zero in the Rashba model (Fig. [3]:). Upon tuning the 
Rashba splitting continously in our model, we find that the lifetime enhancement factor increases 
approximately exponentially with Rashba splitting - a consequence of the exponential behavior of 
carrier occupation numbers near the tail of the Fermi-Dirac distribution. Our model predicts that a 
Rashba spliting of 0.1 eV can give rise to an order of magnitude enhancement of carrier lifetime. 

As we have seen, the presence lifetime enhancement is enabled by the relative spin helicity 
of VBM and CBM, while the amount of enhancement depends on the magnitude of the Rashba 
splitting. The symmetry-breaking distortions that influence relative spin helicity and splitting mag¬ 
nitude are therefore intimately related to the SOC enhancement of carrier lifetime. To reveal the 
relation, we start with a tight binding model of the inorganic Pbl^ lattice (see Methods). The 
displacements of Pb atoms along the z direction give rise to effective hoppings between Pb ,v- and 
p-orbitals along the equatorial direction (Fig. [4^), which would vanish by symmetry in the absence 
of such displacements. Similarly, displacement of the apical I atoms along the z direction changes 
Pb-I bond lengths and effective hoppings along the apical direction (Fig. [4]b). These modifications 
of hopping parameters create the effective inversion symmetry breaking electric fields described 
by Hvr in our low-energy theory model (Eq. 1). In this model, we find that the spin textures of 
the valence bands and conduction bands depend on the combination of Pb and I displacements. 
By shifting Pb and I atoms, we can control the spin textures of the valence and conduction bands, 
creating favorable and unfavorable relative spin helicities. As we now proceed to show, this picture 
is confirmed with DFT calculations. 

We focus on two phases of MAPM 3 , pseudocubic phase (a phase) and tetragonal phase (/3 
phase )]- 11 The pseudocubic phase of MAPM 3 has space group P4mm with a — b — 6.31 and c — 
6.32 , which can only hold ferroelectric distortions, such as displacements of Pb and apical I atoms. 
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Shown in Fig. |4]i are the DFT-calculated spin textures and averaged Rashba splittings at VBM and 
CBM respectively for a given pair of Pb and I displacements (see Methods). The areas outlined 
by the solid red lines indicate the structures with favorable relative spin helicities, which also 
have relatively large band splittings. We find that these structures have Pb and apical I displaced 
along opposite directions. This is consistent with typical ferroelectric distortions in inorganic 
ferroelectrics such as BaTi 03 and PbTi 03 . It is noted that large apical I displacement tends to 
drive the system away from the region with favorable relative spin helicity to unfavorable spin 
helicity (Fig.[4]T). 

We further explore the relationship between relative spin helicity and local distortions in tetrag¬ 
onal MAPbl 3 , n which is observed at room temperature. The space group of the tetragonal phase 
is identified as 74cm, allowing both ferroelectric distortion and Pblg octahedral rotation (Fig. [4je) . 
Fig. |4]f shows the spin textures and the averaged Rashba splitting for different Pb and I displace¬ 
ments in 74cm space group. The tetragonal phase exhibits similar displacement-helicity relation¬ 
ship to the cubic phase, indicating that the Rashba splitting can also enhance the carrier lifetime in 
tetragonal phase. 

At room temperature, switchable, large (^100 nm) ferroelectric domains have recently been 
observed experimentally in MAPM 3. 37 Various theoretical studies also suggested the existence of 
local polar regions at room temperature . 24 26 38 The structures that are energetically accessible at 
room temperature are highlighted by dashed red lines in Fig. |4]d and f, which cover a large re¬ 
gion displaying favorable spin textures. Hence, our proposed mechanism is likely to be realized in 
MAPM 3 at room temperature, and provides a possible explanation of long carrier lifetime. More¬ 
over, recent studies of carrier dyanmics under magnetic field also illustrate the significance of spin 
in carrier recombinations . 39 ' 41 

Rashba SOC enhanced carrier lifetime highlights the potential of 3D Rashba materials for 
photovoltaic applications. The ability to incorporate different organic molecules in OMHPs pro¬ 
vides a robust avenue to design 3D Rashba materials. By tuning the dipole magnitude of organic 
molecules in halide perovskites, the spin helicities and band splittings can be manipulated via the 
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Pb and I displacements as discussed above. Conventional experimental techniques of controlling 
bulk polarization (e.g, epitaxial strain) can also be applied to optimize power conversion efficiency. 

In summary, we have proposed an intrinsic mechanism for enhancing carrier lifetime in 3D 
Rashba materials. In the case of OMHPs, such mechanism can be realized by the joint action of 
molecules (electron-phonon coupling) and PM 3 sublattices (giving rise to spin-orbit coupling). The 
photoexcited carriers quickly relax to band edges due to the electron-phonon coupling. When the 
spin textures for CBM and VBM are opposite, the Rashba splitting of bands close to the band gap 
results in spin-allowed and spin-forbidden recombination paths. The spin-forbidden recombination 
path has slow transition rate due to mismatch of spin and momentum. The spin-allowed recombi¬ 
nation path, though kinetically favorable, will only influence a smaller amount of carriers. In order 
to achieve this favorable spin helicity, we explore different Pb and I displacement giving rise to dif¬ 
ferent spin textures. A tight binding model is developed to explain this spin-displacement relation. 
This mechanism allows OMHPs to behave like direct-gap semiconductors upon photoexcitation, 
and like indirect-gap semiconductors during radiative recombination, simultaneously harnessing 
the large carrier densities of the former and the long lifetimes of the latter. Our work is the first to 
introduce this mechanism and highlights the importance of Rashba effect and structural distortion 
for achieving long carrier lifetime and consequently long diffusion length in organometal halide 
perovskites. 

Methods 

DFT and Electron-phonon coupling The plane-wave DFT package QUANTUM-ESPRESSO ^- 42 
with the Perdew-Burke-Ernzerhof 43 functional is used to perform electronic structure and electron- 
phonon coupling calculation. Norm-conserving, designed nonlocal pseudopotentials were gener¬ 
ated with the OPIUM package . 4445 We calculate the inelastic phonon scattering rate for electrons 
and holes using Fermi’s golden rule : 35 46 


T (ph)-i =h - 
^k.n^n 


V 


(2k)- 


I / d 3 kV)U 2 m, n -E k ' y 


h(o v 


(4) 
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where k and n are the wavevector and band index respectively, v denotes the phonon mode, E^ n 
are electronic band energies and V is the volume of the unit cell. The phonon frequencies co v and 
electron-phonon matrix elements g n y n k k , are obtained from density functional perturbation theory 
(DFPT)P 

Lifetime enhancement factor The rate constant B x x r in Eq. 3 of spin-allowed (x = x! — ±1) 
and spin-forbidden (x 7 ^ X') transitions are obtained from averaged DFT calculated oscillation 
strength. For example, B u ~ 

Tight-binding model The tight-binding model is based on Pblj~ structure. Pb 6,v, Pb 6 p and the I 
5 p orbitals are included in this model with spin degree of freedom. The tight-binding Hamiltonian 
is: 

#tb = #ho P + tfsoc (5) 

where H bop considers the nearest neighbor hopping between two orbitals, described by t sp , t ppa and 
tppn for s-p o hopping, p-p a hopping and p-p n hopping respectively. Hsoc 1S the on-site SOC 
term defined as AR as hb a L • S. Hopping parameters and Apashba are fitted to DFT band structures of 
pseudocubic MAPM 3 with experimental lattice constants. We reduce our tight-binding Hamilto¬ 
nian into the Rashba effective model (Eq. 1) in two steps, following the procedure outlined in . 28 
First, the I p orbitals are removed by projecting Hj b to the subspace of Pb orbitals: 

H Pb = 1 ( 6 ) 

E — £2Hj B £2 

where & and XX. are projection operators to the Pb and I subspaces respectively. This results in 
an effective Hamiltonian H Bb containing inversion symmetry breaking terms which modify the 
effective hopping between Pb orbitals . 28 Next, H Bb is reduced to Hr using a similar projection to 
the CBM and VBM. 

Rashba splitting energy phase diagram The Rashba splitting energy is defined as |A£^| = 
min [<£c,*=cbm+i “^c,*=cbm> , <£v, j= vbm -£v,j=vbm-i)] • Here, the sign of A E R is indicated 
as +1 for favorable spin helicity, and —1 for unfavorable spin helicity. “()” indicates the average 
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over k points near CBM or VBM. 
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Figure 1: Diagram of Rashba bands and the electron transport path. The cyan and orange arrows 
indicate the directions of the spins. The spin texture % indicates spin vortex direction, with its 
signs characterizing spin rotation in “clockwise" (x — — 1) and “counterclockwise" (x = +1). 
After absorbing the photons, the excited electrons on conduction bands C x=+ i and C^=_i will 
quickly relax to C x= -i band minimum due to the inelastic phonon scattering. Similarly, the holes 
will quickly relax to the V^ =+ i band maximum. However, the radiative recombination of C^=_i 
—> V%=+i is a spin-forbidden process due to the opposite spin states they have. Moreover, the 
minimum of C x= -i band and the maximum of Vy =+ i band are located in different positions in 
the Brillouin Zone. This creates an indirect band gap for recombination, which further slows down 
the recombination process. 
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Figure 2: Dependence of phonon-induced relaxation rate on carrier energies (blue lines) for elec¬ 
trons (positive energies) and holes (negative energies). The VBM is located at E — 0 eV, and the 
CBM is located at E = 0.73 eV. The energies of phonon modes that contribute strongly to carrier 
relaxation are shown as dotted lines. 
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Figure 3: a) Population of carriers at T -298 K, calculated from the Fermi-Dirac distribution, for 
a) the top valence band and b) the bottom conduction band of MAPM 3 . Shown in dashed lines 
are the population of carriers in a model material with the same band masses as MAPM 3 , but with 
vanishing Rashba splitting, c) Unitless lifetime enhancement factor, as defined in the text, as a 
function of the conduction band Rashba splitting energy. 
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Figure 4: a) and b) Hopping schemes illustrate the effective electric field Hvr jSP and hvR, pp created 
by vertical hopping and horizontal hopping respectively for a range of Pb and I displacements. 
These two factors caused by different Pb and I displacements controls spin textures of CBM and 
VBM differently, giving rise to different spin helicities. c) Schematic diagram showing Pb and 
I displacement in pseudocubic MAPM 3 . Pb: Silver. I: Indigo. Broken circles are original high- 
symmetry positions. Molecules are not shown here, d) Phase diagram of splitting energy and 
spin texture for structures with different Pb and apical I displacement in pseudocubic MAPM 3 
calculated from DFT. The color is the minimum value between the averaged splitting energy of two 
Rashba conduction bands and two valence bands (see Methods). The spin texture phase boundaries 
are indicated by the solid red lines. When the structure transforms from a favorable spin texture 
region to an unfavorable spin texture region, the two Rashba valence bands or conduction bands 
exchange, creating negative splitting energy. The dashed lines indicate the areas with energy cost 
less than 25 meV (under room temperature fluctuation) to distort Pb and I. e) Schematic diagram 
showing Pb and I displacement in tetragonal MAPM 3 . f) Similar to d, phase diagram of splitting 
energy and spin texture for structures with different Pb and apical I displacement in tetragonal 
MAPM 3 calculated from DFT. 
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